BCS BEC crossover and phase structure of relativistic system: a variational approach 
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We investigate here the BCS BEC crossover in relativistic systems using a variational construct 
for the ground state and the minimization of the thermodynamic potential. This is first studied 
in a four fermion point interaction model and with a BCS type ansatz for the ground state with 
fermion pairs. It is shown that the antiparticle degrees of freedom play an important role in the BCS 
BEC crossover physics, even when the ratio of fermi momentum to the mass of the fermion is small. 
We also consider the phase structure for the case of fermion pairing with imbalanced populations. 
Within the ansatz, thermodynamically stable gapless modes for both fermions and anti fermions are 
seen for strong coupling in the BEC regime. We further investigate the effect of fluctuations of the 
condensate field by treating it as a dynamical field and generalize the BCS ansatz to include quanta 
of the condensate field also in a boson fermion model with quartic self interaction of the condensate 
field. It is seen that the critical temperature decreases with inclusion of fluctuations. 

PACS numbers: 12.38.Mh, 24.85,+p 



I. INTRODUCTION 



Color superconductivity has become an active field of 
research during the last few years in the field of strong in- 
teraction physics [1] . At asymptotically high density the 
ground state of QCD is shown to be a color supercon- 
ductor from first principle calculations [2]. However, for 
intermediate densities, relevant for the matter in the core 
of neutron stars, the perturbation theory breaks down 
and one uses effective models of strong interactions whose 
parameters are fitted to reproduce low energy hadronic 
properties. Although, it is reassuring that certain quan- 
tities like the superconducting gap estimated from weak 
coupling QCD or using effective models like Nambu Jona 
Lasinio (NJL) model yield similar magnitudes, the strong 
coupling dynamics at not so high densities remain un- 
controlled. As the density decreases, before the quarks 
are confined, the coupling can become large enough so 



* Electronic address: bhaswar@prl.res.in 
i Electronic address: hm@prl.res.in 

* Electronic address: amruta@physics.iitd.ac. in 



that the coherence length in the superconducting phase 
can be of the order of inter particle separation [3, 4]. 
When this happens it is natural to imagine the quark 
quark Cooper pairs as localized bound states rather than 
describing them as an extended macroscopic medium. 
At low enough temperature the ground state can then 
become a Bose Einstein condensate (BEC) of diquark 
molecules. At still lower densities this diquark matter 
can undergo a confinement transition to form hadronic 
matter. The scenario that we have in mind is similar 
to as in Ref. [5], where, as the baryon density is in- 
creased, the baryonic matter undergoes a phase transi- 
tion to diquark BEC phase. The transition from BCS 
to BEC is most likely to be a crossover similar to their 
non relativistic counterparts in condensed matter sys- 
tems like cold fermionic atoms. Various methods have 
been used recently to describe relativistic BCS - BEC 
crossover within a four fermion point interaction model 
[4, 6, 7]. Effects of fluctuations have also been studied 
in the recent past to go beyond the mean field approxi- 
mations [8-11] using a gaussian approximation as well as 
including collective excitations [12]. Another interesting 
feature for the matter in the interior of neutron star is the 
variety of exotic non BCS phases that arise when kine- 



matical constraints like neutrality with respect to color 
and electric charges are imposed. Such kind of stressed 
pairing has attracted attention both in quark matter [13] 
and in polarized Fermi gas of atoms [14, 15]. Relativistic 
BCS-BEC crossover was also studied in a boson fermion 
model and a rich phase structure was observed in such a 
system when there is a mismatch in the number densities 
of the condensing fcrmions [6] . 

Our approach to such problems in quark matter [16- 
18] as well as in cold atoms [19] has been variational. 
Charge neutral matter was considered with an explicit 
construct for the ground state in terms of quark pair op- 
erators as well as quark antiquark operators. The ansatz 
functions in this construct were determined from the min- 
imization of the thermodynamic potential subject to the 
constraints of neutrality conditions with respect to color 
and electrical charges. The effect of six quark determi- 
nant interaction was also considered which highlighted 
the consequences of strange quark mass in such stressed 
pairing cases [20] . We apply here similar method to study 
BCS-BEC crossover (as well as possible phase transition 
for mismatch in number densities) for a relativistic sys- 
tem in a fermionic theory with a four fermi point inter- 
action model. As we shall demonstrate here the simple 
BCS type of ansatz for the ground state leads to results 
similar to the mean field results of a relativistic Bose 
Fermi model [6]. This is demonstrated by taking a BCS 
type ansatz for the ground state and treating the con- 
densate field as a classical auxiliary field. We next treat 
the condensate scalar field as a dynamical field and gen- 
eralize the BCS ansatz to include the quanta of this con- 
densate fields in a Bose Fermi model with a quartic self 
interaction for the scalar field. This improved ansatz for 
the ground state leads to a mass gap equation for the 
scalar field and the corresponding thermodynamic po- 
tential that can be obtained by resummation of bubble 
diagrams of perturbation theory similar to the results 
obtained in Cornwall- Jackiew Tombolis (CJT) compos- 
ite operator formalism [21, 22]. 

We organize the paper as follows. In section II we 
shall consider a relativistic model with a four Fermi point 
interaction to discuss the BCS-BEC crossover physics. 
We shall take the bosonic condensate field as a classi- 
cal auxiliary field while retaining the quantum nature of 
the fermionic field. In subsection II A, we spell out the 
ansatz for the ground state and in subsection II B, we 
evaluate the thermodynamic potential and consider the 
minimization of the same with respect to the ansatz func- 
tions. In section III, we solve the gap equation and the 
number density equations and evaluate the thermody- 
namic potential. The thermodynamic potential for dif- 
ferent phases are compared in this section to decide which 
phase is thermodynamically stable at what coupling and 
at what difference in the chemical potential. We also dis- 
cuss our numerical results in this section. The results ob- 
tained in this section using the variational ansatz, is simi- 



lar to that obtained within mean field approximations in 
earlier calculations regarding BCS -BEC crossover and 
phase structure for imbalanced populations [4, 6, 7]. In 
section IV, we treat the scalar condensate field as a dy- 
namical field and generalise the ansatz of section III to 
include the quanta of the fluctuations over and above the 
BCS ground state. We investigate this in a model with 
quartic self interaction for the scalar field, calculating the 
mass of the scalar field self consistently. We also include 
the effect of temperature in these calculations. Finally 
we summarize our results in section V. 



II. A FERMIONIC MODEL FOR 
RELATIVISTIC SUPERFLUID 

We shall consider a general relativistic model with only 
fermionic degrees of freedom in the Lagrangian. In par- 
ticular we consider a general form for the Lagrangian 
given as 

C = f( l7 M 9, - m + w a W + C I = C } + C I (1) 

where, ip l ,tp l denote the Dirac fields with indices i = 
1,2 denoting the 'flavor' for the fermions. Further, we 
have taken, for simplicity, both the fermions with the 
same mass m. With different masses, the fermi energies 
or the chemical potential in the weak coupling, can be 
different for the same number densities of the two pairing 
species. In the present analysis, for simplicity, we assume 
the masses of the pairing fermions to be the same, but 
consider their chemical potentials to be different. 

To study superfluidity that results from Cooper pair- 
ing of two different flavors of fermions, we introduce the 
following local interaction term in the Lagrangian, 



Ci 



-G^V'l^'lX^SVcl^l 



(2) 



Here, G is the coupling constant, ip c = C-tp T , V> c = ip T C, 
with the charge conjugation matrix being C = ij 2 j°. 
The |e 4J | term ensures the cross flavor, spin zero anti- 
symmetric pairing. 

To illustrate the mean field calculation, introducing a 
field $, we can rewrite the interaction term as 



£/ = g\e tj \ (V'VV^ + V'cTV** 



m^<r$ (3) 



with G — g 2 /m 2 which can be identified upon elimina- 
tion of the field $. We note that in the absence of a 
kinetic term the field $ is an auxiliary field whose value 
in terms of the fermion bilinears is known. In the fol- 
lowing, however, we shall treat the field $ as a classical 
field and replace it by its expectation value while retain- 
ing the quantum nature for the fermion field. This will 
enable us to calculate the effective potential as a function 
of 0o = ($>■ 



A. The ansatz for the ground state 

To make the notations clear, let us first note that the 
fermion field operator expansion in momentum space is 
given as [24, 25] 



V(x) 
1 

W) 



1 



V>(k)e 2k ' x dk 



(2tt) 3 / 2 
j-j^ J [J7 (k)g(k) + V {-k)q{-k)) e^dk (4) 



where, 





-a ■ fcsin(— 
cos(4j- 




(5) 



The operators q and q are the two component particle an- 
nihilation and antiparticle creation operators respectively 
which annihilate or create quanta acting upon the per- 
turbative vacuum |0). We have suppressed here the flavor 
indices of the fermion field operators. The function x°(k) 
in the spinors in Eq.(5) are given as cotx? = m.j/|k|, for 
free massive fermion fields, i being the flavor index. For 
massless fields % (|k|) = ir/2. We shall in the following, 
however consider for simplicity, toi = nri2 = m. We shall 
use the notations and conventions of Ref. [17, 18] and 
recapitulate briefly the construction of the variational 
ansatz for the ground state. We take it as a squeezed 
coherent state fermion condensates given as [16-18] 



|fi)=W«i|0). 



(6) 



Here, Lid is a unitary operator which creates fermion 
pairs. Explicitly, the operator, hid is given as 

U d = cxp(^ - B d ), (7) 

where, B d is the pair creation operator as given by 



Bl 



«*(k)tr/(k)gi r (-k)t|< 
'#(k)r/i(k)e r (-k)h 



dk 



dk. 



(8) 



In the above, i,j are flavor indices, and r(= ±1/2) is the 
spin index. Further, /(k) and /i(k) are ansatz functions 
associated with fermion pairs and anti fermion pairs de- 
scribing the condensates and shall be determined by the 
minimization of the thermodynamic potential. Note that 
we have assumed these 'condensate functions' /(k) and 
/i (k) to be independent of flavor indices. We give a post- 
facto justification for this to be that the functions depend 
upon the average energy and average chemical potential 
of the fcrmions/anti-fcrmions that condense. 



Finally, to include the effects of temperature and den- 
sity we write down the state at finite temperature and 
density \il(/3, /J,)} taking a thermal Bogoliubov transfor- 
mation over the state |f2) using thermo field dynamics 
(TFD) as described in Ref.s [26, 27]. We then have, 



where Up^ is given as 



M^UdlO), 



U 



/3,M 



B\p,v)-Btf,v) 



(9) 



(10) 



with 



B<<(/3,n) 



g , (k)t^_(k,/3,A*)«'(k)t]dk 

[q'(k)6 + (k,(3,n)q;(k)\dk. (11 



In Eq.(ll) the ansatz functions 0±(k, (3, /z) will be related 
to quark and antiquary distributions, and, the under- 
lined operators are the operators in the extended Hilbcrt 
space associated with thermal doubling in TFD method. 
In Eq.(ll) we have suppressed the flavor indices in the 
fermion operators as well as in the functions #(k, /3,/x). 

All the functions in the ansatz in Eq.(9) are to be ob- 
tained by minimizing the thermodynamic potential. We 
shall carry out this minimization in the next subsection. 



B. Evaluation of thermodynamic potential and gap 
equations 

To calculate the thermodynamic potential correspond- 
ing to the Lagrangian of Eq.(l) and Eq.(3) and the state 
given in Eq.(9), we first write down the expectation val- 
ues of the following fermion operator bilinears. 



and, 



(n(p,n)\^(k)^ s (ky\n(p,n)) 

= ^AV 7 ,(k,/3, M )<S(k-k') 



(n(/?,,u)|$ + (k)$(k')|fi(M> 

= S^A% 5 (k,p^)S(k-k') 



(12) 



(13) 



where, 



A* ±7 ,(k, 0,vi) = -[l± (F*(k) - F*(k)) ± ( 7 ° cos X 4 (k) 



a ■ ksinx^k)) (l - F l {k) - F^k)) 



7 5 



(14) 



In the above, ip(k) is the Fourier transform of ^(x) [18]. 
The effect of the fermion condensates and their tempera- 
ture and/or density dependences are encoded in the func- 
tions F l (k) and ^(k) given as 



(15) 



F l (k) = sin^ 9 l _ (k) + &\tl z f (k) cos 28 1 ! 1 (k) 



and, 



F^k) = sin 2 6\ (k) + sin 2 /i (k) cos 2d% 3 (k) . (16) 



We have defined cos 20± = 1 — sin 2 9 l ± — sin 2 6 3 ± with 

i + 3- 

For difermion operators, we have, 



M/MIC(x)v4(°)in(M> 



Ue^JM,^ 



and, 



(27T) 



(QtMIV^W^WIfiCM) 



(2tt) 3 



ikx-r,i,j 



Pi' (k,/3,//)dk, 



where, 



(17) 



(18) 



7^(k,/3, M ) = V 



S' ( k)(co S (^i _ 7 .k sin *L^ 



^(k,/3, M ) 



o / X» + Xj 
7 cos - 

1 \ 2 



|e y '|C 7 5 



-a-kshi [ Xi + Xi ) U«(k 



+ 7 U cos 



5 4J (k) (cos 
X* + Xj 



Xi Xj 



7 • k sin 



1 5 C 

Xi ~ Xj 



a-ksh.l Xi + XJ ) U«(k) 



r 



(19) 



(20) 



Here, C = «7 2 7° is the charge conjugation matrix (we where e is the energy density and s is the entropy density 

use the notation of Bjorken and Drell) and the functions and p % — (ip^tj; 1 ) (i = 1, 2) is the number density of i-th 

5(k) and A(k) are given as, species. 

S ' W = sin2/(k)cos20_l (k, 0, n) j^ ; s ^hen straightforward to calculate the expecta- 

+ sin2 fi ( k) cos 20 1,3 ( k u) (21) tion va l ue of the Hamiltonian corresponding to the La- 

grangian given in Eq.(l) and Eq.(3). This can be written 

and, as 



A^{k) = sin2./ , (k)cos2(f J (k,/3,^) 

- sin2/ 1 (k)cos2^'(k,/3, M ), (22) 

These expressions are used to calculate thermal expecta- 
tion value of the Hamiltonian and compute the thermo- 
dynamic potential given as 



fi = e - mV - -j s 



e - M y = (H- M W) =T + V D 



(24) 



Explicitly, the kinetic energy minus the /i' p l part is given 
(23) as 

I 



2 

T = (Sl(J3,n)\il>U-i<* • V + 7 V - m'M W,A«)) = -^ ^ J dk (Vk 2 +m 2 (F* + F{) - p\F* - F{ 



(25) 



where, F l and F[ are given by equations (15) and (16). 
Here we have subtracted out the vacuum contributions. 



Similarly, the contribution from the interaction from 



Eq.(3) to the energy density is given as 
V D = -(n{0,n)\C I \iW3,n))=-4gI D <h+rr%<ft (26) 



where we have taken <bn to be real. In the above, 



J 



In = ^cY\e iJ W)- {2 



i-3 f dk [sin2/(k)(l - sin 2 9 l _ - sin 2 9 2 _) + sin2/i(k)(l - sin 2 Q\ - sin 2 $%)] 



r 



(27) 



which is proportional to the fermion condensate. 

Finally to calculate the thermodynamic potential we 



have to include the entropy density for the fermions. This 
is given as [26] 



J 



(2 



~^E /dk(sm 2 (9imsii 



2 ai 



-,2 ai 



cos 2 9'_ In cos 2 9 l _ + sin 2 6\ In sin 



1 ai 



,2/3i 



cos 6\ In cos 9 



2 ai 



(28) 



r 



The cxtrcmisation of the thermodynamic potential Eq. 
(23) with respect to the condensate functions /(k) and 
/i(k) yields 



tan2/(k) 



2ff0o 



and 



tan2/i(k) = 



2g0o 
e + /U 



/i 



(29) 



(30) 



where, we have defined the superconducting gap A = 
2gcf)Q. In the above e = (e± + e%)/2 and p, = {ji\ + H2)/2 
with ej = -\/k 2 + m 2 . It is thus seen that the conden- 
sate functions depend upon the average energy and the 
average chemical potential of the fermions/anti-fermions 
that condense. 

Finally, the minimization of the thermodynamic po- 
tential with respect to the thermal functions 0± (k) gives 



where, 



and 



:(k) 



,(!) 



1 



exp(/3w 



± ) 



w± = u>± + S £ ± (5 P 



(2) - c r 



(31) 

(32) 
(33) 



Here, iu ± = ^ + ^1, £± = (£i± + 6±)A &± 
ei± fii. Further, the chemical potential difference 5 M 



(/■H - A*2)/2 and 5 e = (ei - e 2 )/2 

and for equal masses for the two species 8 f — with 



Notice that for 5^ > 



m 

(i) 



we can have the possibility of 

(2) 



ei = e 2 = e = Vk 2 
gaplcss modes for wr' or u> 

Using these dispersion relations, and substituting the 
condensate functions and the distribution functions, 
leads to the thermodynamics potential given by equation 
(23), as 



n 



(2tt) 3 
2 

(2tt) 3 /? 

ln(l + 



(2e 
E 



idk 



dk 



In 1 



,(-0"-) 



-M°) 



m 2 ^ 2 ,. 



(34) 



Here the extremisation over </> is yet to be done. An ex- 
trcmisation with respect to <j>q leads to the gap equation 



4 3 2 



dk 

(2tt)' 



cos 29 _ 



W- 



cos 26>_ t 

U1+ 



(35) 



with cos 26>^' 2 = 1 - sin 2 6^ - sin 2 Q\. 

The gap equation is quadratically divergent which is 
rendered finite in the NJL model with a momentum cut 
off A. In the non relativistic case this is rendered finite by 
subtracting out the vacuum contribution and relating the 
four fermion coupling to the s-wave scattering length [19, 
28]. A similar approach can be done for the relativistic 
case also by relating the coupling to the s-wave scattering 
length. This leads to the renormalized gap equation [7] 



m 

4ira 



dk 

(2^)~ 3 _ 



cos 26 r 



u;_ 



cos 29^ 



m 



(36) 



However, after this subtraction, unlike the nonrelativis- 
tic case, the ultraviolet cutoff dependence is still there 
in the above gap equation, although the dependence 
is milder. We might also note that one could have 
defined a renormalizcd boson mass mb, r with m 2 r = 
dfl/d(f>Q\,j >a= T=o,fj,=m as in Ref. [6] and one would have 
arrived at the same gap equation. For the present, we 
shall take the renormalized coupling as in Eq.(36) and 
treat this as the crossover parameter. As a function of 
this coupling one has to calculate the gap parameter A 
for different densities of the fcrmions of the two species. 
The average number density is given as 



P = 



Pi +P2 



= P- ~ P+, 



where, the fermionic component is given as 

i r e- 



p- 



(2^) 3 J u 



cos 20 ' dk 



(37) 



(38) 



and the anti fermionic component of the average number 
density is 



1 



(2tt) 3 7 LU+ 



^ + cos29]i 2 dk 



(39) 



where, cos26^' 2 = (1 - sin 2 6^(k) - sin 2 0|(k)) with 
sin 9±(k) being the thermal distribution functions for 



the fcrmions defined in Eq.(31). The difference in the 
number densities is given as 



Pi - Pi 



1 



[(sin 2 9l- sin 2 e\) 



2 (2tt) 3 

- (sin 2 9 2 _ - sin 2 9%)] dk 



(40) 



p l = [ip 1 ip 1 ) , where the expectation value is taken with 
respect to the state given in Eq.(9). 

Using the gap equation Eq.(35), the thermodynamic 
potential given by Eq.(34) can be rewritten as 



H(A,m,^,/3) 

2 



rfk 



£--«- 



A 2 
2u>_ 



■€,+ ~U-\ 



(2tt) 3 /3 
■ 2 { In (l + e<-*"- >) + In (l + e^^) 



2w 4 



(41) 



To compare the stability of various phases we compare 
the thermodynamic potentials of these phases with re- 
spect to that of normal matter. This can be obtained 
from Eq.(41) in the limit of gap A — > 0. We consider 
the difference in the thermodynamic potentials between 
condensed phase and the normal matter as given by 



0(A, fi, S„, /3) = fi(A, ji, 6^,(3)- fi(A = 0, fi, &„, (3) 
2 /•/,,, A 2 



(2tt) 3 

2 



IM-«- 



(27r) 3 /3 ^ 



al>2 



2u>- 



cos20iT + |£+| -w H 



A ■cos20i ,2 '| dk 



2w 



^ / {ln(l + e (-^-))+ln(l + e (-^+))}-{ln(l + e(-^o-))-ln(l + e ^<+))} rfk(42) 



r 



In the above, Wq-^ = |£ T | =F 5, 



w, 



(2) 

0=F 



|£fI±<^> 



spond to the normal matter dispersion relations for the 
two species. For stability of the condensed phase f2 has to 
be negative with A and p, determined from the gap equa- 
tion Eq.(36) and the number density equation Eq.(37). 
Further, one has to ensure that the solution corresponds 
to a minimum and not a maximum. In what follows we 
shall restrict ourselves to the case of zero temperature 
only. As noted earlier we shall consider here without 
loss of generality, <5 M > 0. Also we consider a system 



of equal masses for the fermions so that S e = 0. This 
leads to the possibility of quasi particle energy for species 
'1', u_ or the quasi anti-fermion energy for species '2', 

(2) 

w\ becoming negative. In that case the distribution 
functions given by Eq.(31), become Heaviside (O) - func- 
tions i.e. sin 2 9 a = 0(—ui a ). Further, using the identity 
linia^oo ln(l + cxp(— ax)) / a = —x<3(—x), in Eq. (42), 
the zero temperature thermodynamic potential becomes 



no(A,jit, jp 



(2tt) 3 



W (1) - 



(2tt) 3 

A 2 



2w_ 



IC- 



(-*a+ 



A 2 
2d>_ 



,(2) 



_A2 

2(1-4 

r 



2(I>, 



,(2) 



dk 



,(1) 



Uh 



,( 2 ); 



,(2)- 



-<') - U^0(-C^) - <'*(-<') 



dk 



(43) 



The equations (36) and (37) need to be solved self 
consistently to determine the gap as a function of the 
coupling and its stability will be decided by calculating 
the thermodynamic potential. We note here that BEC 
is usually discussed using the canonical ensemble where 
the particle number density is fixed as an external pa- 
rameter. We shall also consider fixed number density 
here to discuss BCS-BEC crossover/phase transition in 
the relativistic fermionic system [6, 8-11]. However, we 
might also note here that, to discuss quark matter, one 
usually employs grand canonical ensemble with a fixed 
quark chemical potential to explore the QCD phase di- 
agram in the chemical potential and temperature plane. 
In the numerical calculation that follows we keep the av- 
erage number density fixed and consider the solutions as 
a function of the coupling and the difference in chemical 
potentials. Sometimes we find multiple solutions for the 
gap and average chemical potential satisfying Eq.s (36) 
and (37) corresponding to multiple cxtrema of the ther- 
modynamic potential. In such cases, the solution which 
has the least thermodynamic potential is chosen. We 
also verify for the positivity for the second derivative of 
the thermodynamic potential for this solution. This way 
we ensure that the pair of solution for supcrfluid gap and 
the chemical potential correspond to the minimum of the 
thermodynamic potential. The detailed numerical calcu- 
lations of the present investigation are discussed in the 
next section. 



III. NUMERICAL SOLUTION OF THE GAP 
EQUATION AND PHASE STRUCTURE 

For numerical calculations it is convenient to introduce 
the dimcnsionlcss quantities in terms of Fermi momen- 
tum kf or Fermi energy e/ = 



k 2 , + to 2 , defined through 

|k| = kfX, r\ = l/(fc/a), to = fc/TO , A = efz, fi = efji. 
The gap equation at zero temperature can then be writ- 
ten as 



inn ./o 



dxx 



1 



1 



2e(x) 



^(-^Vt^-^) 



U) 



(i) 



CO 



(2) 



(44) 



Similarly, the equation for the average number density 
given by Eq.(37) can be rewritten in terms of these di- 
mcnsionlcss quantities as 



pxmax 

1.5 / dxx 

Jo 



u>+(x) 



e(-4 2) : 



u>—(x) 



l-Q(-u 



(1)' 



for the average 

i(x) ± aV(T 



density. 
to 2 ), e(x) 



Here, 



= Va 



(45) 

i±(x) = 

l 2 , L0± - 



£ 2 _ + z 2 (l + m 2 ) and finally x max = A/fcj is the 
upper cut-off for momentum in units of Fermi momen- 
tum Kf. Here fi is the average of the chemical potentials 
of the two species in units of Fermi energy. 

To analyze the crossover, let us first consider the sym- 
metric case - namely when the chemical potential dif- 
ference between the two species is zero. In this limit, 



©(- 



,W> 



become zero in both the equations (44) and 



(45). We note that we have essentially three dimensional 
quantities in the problem: the cutoff A, the mass of the 
fermion in and the scattering length a. We further note 
that the dimensional coupling G is bounded above with a 
critical value G C A 2 > 2-7T 2 , beyond which the zero density 
vacuum itself is unstable to form fermion pairs leading to 
a Majarona mass for the fermions. 

We might note here that at zero density, the minimum 
excitation energy for the fermion is its mass to. For a nor- 
mal matter, with a finite chemical potential it is (in — /j). 
In the BEC state, the decay mode of the bound state 
is that of the bound state going to two fermions. The 
threshold for this energy is thus 2(m — /i). The bosonic 
bound state should therefore be stable if this threshold 
energy is positive which in turn means that to > /j. This 
we shall take as our working definition for BEC phase to 
distinguish between the BCS and BEC phase as we in- 
crease the dimcnsionlcss coupling parameter r\ from weak 
coupling BCS (large negative 77) to strong coupling BEC 
phase (large positive n) through unitary limit (77=0). 

For numerical calculations, we choose the parameters 
Xmax — A/fc/ = 50, to = m/kf — 5 and study the results 
by varying the dimcnsionlcss coupling n from negative 
to positive values. The resulting gap parameter and the 
chemical potential are plotted in figures 1 and 2. To show 
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FIG. 1: Gap parameter in units of Fermi energy is plotted as 
a function of the dirnensionless coupling. The dotted line cor- 
responds to the case where anti particle contributions are not 
included. The solid line corresponds to the case with inclusion 
of the antiparticle contributions. In this plot, we have chosen 
A/kf = 50 and m/kf = 5 




FIG. 2: The scaled chemical potential (fi — m)/(e/ — m) is 
plotted as a function of the dirnensionless coupling, r\. The 
dotted line corresponds to the case where anti particle contri- 
butions are not included. The solid line corresponds to the 
case with inclusion of the antiparticle contributions. In this 
plot, we have A/kf = 50 and m/kf — 5 



the contribution of the antiparticle degrees of freedom we 
have plotted the results obtained by solving the coupled 
gap equation Eq.(44) and Eq.(45), without and with the 
antiparticle contributions. 

One might naively expect that for this non relativistic 
kf/m — 0.2 case, the antiparticle channel is suppressed. 
But as may be seen from Fig.l, while such an expecta- 
tion may be the situation for weak coupling BCS regime, 
these contributions become increasingly important as the 
coupling increases. As the coupling increases, the chem- 
ical potential Jx = ji — m decreases and changes sign at 
coupling t] rs 1.04 signaling the BEC regime. To appreci- 
ate the relativistic effects we also consider the case with 
m/kf = 0.67 and A/kf = 3.3 [6] and the resulting gap 
and the chemical potential are shown in Fig. 3. 




FIG. 3: Superfluid gap (solid line) and chemical potential 
(dotted line) in units of Fermi energy as functions of the di- 
rnensionless coupling, r\. 

As we might observe, in the weak coupling BCS limit 
(large negative r\), the chemical potential is given by the 
Fermi energy. With increase in the coupling, the chemical 
potential decreases and becomes negligible as compared 
to the Fermi energy in deep BEC regime of large positive 
value of r\. The gap starts with exponentially small val- 
ues in the weak coupling regime as expected from BCS 
theory and rises monotonically as the coupling increases. 
It becomes of the order of the Fermi energy at around 
the unitary regime 77=0. For 77 = 0, the superfluid gap 
and the chemical potential turn out to be A = 0.3e/ and 
[i = 0.78e/ respectively. We might note here that the 
same turns out to be about /u = 0.37e/ in Ref. [6]. The 
discrepancy between the two can easily be understood 




FIG. 4: Number densities of fermions, p- (solid line) and 
anti-fermions, p+ (dotted line) in units ofk^/Zit 2 as a func- 
tion of dimensionless coupling n. 



by comparing the gap equation in the unitary limit in 
the two cases. The difference lies in the way the renor- 
malization of the gap equation is being done in these 
cases. To compare with the non relativistic results, we 
subtract out the mass terms from both the chemical po- 
tential as well as the Fermi energy. The resulting ratio 
then becomes fl/if = (jj, — m)/(e — m) rs 0.5. In the non 
relativistic fermionic models this value turns to be 0.4 to 
0.5 [15, 29-31]. As r\ increases, at about r\ = 1.68, the 
chemical potential becomes smaller than the mass of the 
fermion and the system goes to the BEC regime. 

As the coupling becomes close to the unitary regime, 
the anti fermion contributions become important. In Fig. 
4, we show the number densities of the fermions p_ and 
the anti-fermions p+ as defined in Eq.(38) and Eq.(39) re- 
spectively. While in the BCS regime the antiparticle con- 
tribution to the number density turns out to be negligible 
as compared to the particle contributions, it becomes in- 
creasingly important near the unitary regime. At very 
large values of rj, the chemical potential becomes negligi- 
ble, the contributions of the particle and the antiparticle 
to the number densities become large as compared to the 
density p and their difference produces a conserved net 
density [7]. 

Next, we consider the case of superfluid with a mis- 
match in the chemical potentials i.e. 8^ ^ 0. In the 
numerical calculations, we keep the average density fixed 
and calculate the average chemical potential and the su- 
perfluid gap using Eq. (37) and Eq. (36) respectively. The 



stability of the solution is analysed by calculating the 
thermodynamic potential. Sometimes, particularly near 
the phase transition region, there are multiple solutions 
of the gap and number density equations. Of these, we 
choose the solution which has the least thermodynamic 
potential. We also numerically verify the positivity of 
second derivative of the thermodynamic potential at this 
pair of A and p.. This way we ensure that the solution 
we get is a minimum and not a maximum of the thermo- 
dynamic potential. 

To analyze the nature of the solutions, we shall fur- 
ther assume, without loss of generality, 5 U > 0. In that 

case, the quasi particle energy for species '1', u_ (k) = 
cD_(k) — Sfj,, and the same for quasi antiparticle energy 

of species '2', w+ (k) = w+(k) — S^ can become neg- 
ative. In that case, at zero temperature the contri- 
butions of the O - functions both in the gap equation 
Eq.(35), and, the thermodynamic potential Eq.(34) can 
be non vanishing. The - functions limit the range 
of the momentum integrations. Solving for the zeros of 
ui]_ , we note that this vanishes at momenta k 2 . 



min max 



(2) 



Similarly ui\ vanishes at mo- 



(a±^-a>) - 

menta satisfying k 2 , , . . _ , 

Clearly this is possible provided the gap A is smaller 
than 5 p. The zeros of the dispersion relations correspond 
to effective Fermi surfaces. In general there can be two 
fermi surfaces for species '1' along with the gapped ones. 
From the dispersion relations for the quasi particles and 
antiparticles, it is clear that we also can have the inter- 
esting possibility of interior gap solutions for particles of 
species '1' and antiparticles of species '2' [6]. 

After these general remarks regarding the topology of 
Fermi surfaces, let us discuss some numerical results of 
the present investigation. As mentioned earlier, we keep 
the average density fixed and calculate the superfluid gap 
and the average chemical potential for a given coupling 
and a given chemical potential difference. In case there 
are multiple nontrivial solutions for the gap and the av- 
erage chemical potential, the solution for the gap and av- 
erage chemical potential pair is accepted which has the 
least thermodynamic potential. 

In the numerical calculations that we present here, we 
take the parameters A/kf = 3.3 and m/kf — 0.67. In 
Fig. 5, we plot the quantity <5„/A, the ratio of maximum 
chemical potential difference to the superfluid gap, that 
can sustain pairing, as a function of the dimensionless 
coupling r\. For a larger value of <5 M , beyond 8 C there 
are no acceptable solutions of the gap and number den- 
sity equation with a non vanishing gap and with a lower 
thermodynamic potential. For weak coupling in the BCS 
limit, it approaches the Clogston - Chandrasekhar limit 
of <$£/ A ~ 0.72. Initially, as the coupling increases from 
BCS to BEC regime, 6°/ A increases monotonically as 
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FIG. 5: Ratio of critical chemical potential difference to the 
gap as a function of the coupling strength rj. Gap less phase 
appear for r\ > 1.9. Solid line denotes the BCS regime and the 
dotted line indicates the regime where quasi particles of species 
'1 ' become gapless. The dashed line indicates the regime where 
the antiparticles of species '2' also become gapless. 



FIG. 7: Dispersion relation for gapless modes. Solid line 
shows the dispersion relation for quasi particle of species 
'l'((jj_ (k)) and the dotted line shows the dispersion relation 
for quasi- antiparticle of species '2' (u>\ (k.)). The plot is for 
the case of coupling r\ = 3 and <5„/Aq = 1.195. 




i.O 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
k/k F 



FIG. 6: Quasi particle dispersion relation for species '1 '. The 
plot is for the case of coupling r) = 2.1 and chemical potential 
difference S/Ao = 1.125, Ao being the gap at zero chemical 
potential difference. 



shown by the solid line in figure 5. In this regime there 
are no gapless modes and the density difference, S p , be- 
tween the two species is zero. At r\ ~ 1.9, 5 C / A increases 
sharply with rj until it reaches a value of the order of 2.2. 
This is shown by the dotted line in Fig. 5. In this regime, 
lu_ becomes gapless while all other modes are gapped. 
In this region, (p, — m) is negative and hence this gapless 
phase lies in the BEC regime. The density difference be- 
tween the two species becomes nonzero. The dispersion 
relation ui_ (k) in this region is shown in Fig. 6. In par- 
ticular we have chosen here, i] = 2.2 and S^/Aq — 0.47. 
Here Ao is the superfluid gap at zero chemical potential 
difference and turns out to be A = 0.637e/. With these 
parameters, the average chemical potential turns out to 
be Jl = 0.47e/ and the gap is A ~ 0.35e/. 

The ratio S c /A increases with rj till rj ~ 2.5. Beyond 
this, there are two gapless modes. Both the usual quasi 
particle for species '1' and the quasi anti particles of 
species '2' become gapless beyond this coupling. The 
dashed line in Fig. 5 represents this regime. In Fig. 7 we 
show the dispersion relation for these gapless modes. We 
have taken the coupling r\ = 3 and 5/Aq — 1.195 here. 
The values for the gap and the average chemical potential 
here is A = 0.72e/ and p, = 0.23e/ respectively. 

We also would like to note here that we did not find any 
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FIG. 8: Superflmd gap as a function of difference in number 
densities of the condensing species. This is plotted for r\ = 2.1. 
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FIG. 9: Number density distribution for the two species. Solid 
line corresponds to species '1 ' and the dashed line corresponds 
to species '2'. This plot corresponds to r\ = 3.0 and 5 M /Ao = 
1.195 



breached pairing, i.e., two Fermi surfaces with nonzero 
value of the gap for any value of the coupling rj. These 
results are similar qualitatively to those of Ref. [6] where 
a different renormalization was adopted for the crossover 
parameter. 

We might mention here that we do not keep the den- 
sity difference fixed in our calculations. We perform the 
calculations for a fixed average density and a given chemi- 
cal potential difference. In presence of gapless phases, the 
density difference becomes nonzero. In Fig. 8, the de- 
pendence of the gap on the density difference, 5 p is shown 
for coupling rj = 2.1. Superconductivity is supported for 
a maximum density difference of 5 P ~ 0.9(fc?/37r 2 ) be- 
yond which the system goes over to normal matter with 
zero gap. For coupling rj < 1.9 we do not find any su- 
perfluid phase free energetically favorable with any non 
zero value of density difference S p . However, a chemi- 
cal potential difference can still support a Cooper paired 
state. 

Next, let us consider the number density distribution 
of the two species in the momentum space when gapless 
modes exist. Let us note that the number densities at 
zero temperature for the two species are given by 



Pi 



^-3 J dkf sin 2 /(k) + e(-c W ) cos 2 /(k) 



(2 



(2) 



~ (2tt) 3 
for species '1' and 



sinVi(k)(l-^(-^)) 



dkni (k) 



(46) 



P2 



(4i>2 

2 

J2nf 



Vh 



dk[ sin 2 /(k)(l-9(-w w )) 



(sin 2 /i(k)+cos 2 /i(k)0(-4 2) )) 



dkri2(k.) 



(2tt) 3 



(47) 



for species '2'. Here, tan2/(k) = A/<£_ and tan2/i(k) = 
A / £+. In Fig. 9 we have plotted the momentum 
space density distributions ni(k) and ri2(k) for 77 = 3.0 



,(!» 



and 

(i) 



.(2) 



and <5„/Ao = 1.195. In this case both 

become gapless. In the region where both ui y _^' and 



,(2) 



are negative, n\{k) = 1 and ri2(k) = -1 as may 
be seen from equations (46) and (47). In the region, 
where only u>_ is negative, n\{k) = (1 + £ + /tI> + )/2 and 



(1 + £ + /«+)/2. Finally, when both u>_ and 

are positive, both ni(k) and ri2(k) are identical and 
are given by the relativistic BCS distribution function 



n 2 (k) 

,,( 2 ) 
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n(k) = (£ + /w + — £_/u)_). This is precisely what is de- 
picted in Fig. 9. Let us note here that although the in- 
dividual distribution functions in the momentum space 
for the two species could be negative, the average oc- 
cupation number densities n(k) = (m(k) + «2(k))/2 
as well as the difference in occupation number densities 
8 n (k) = (rti(k) — ri2(k))/2 are always positive definite. 

Most of the results obtained in this section for the 
model with a four fermion interaction are similar to the 
mean field results obtained in a boson- fermion model [6] . 
It is nice to see the similarity to the mean field results of 
Ref. [6] which in our investigation arises with a simple 
ansatz for the ground state given by of Eq. (9) deter- 
mined through an extrcmisation of the thermodynamic 
potential. As emphasized in section II, the scalar con- 
densate field was considered as a classical auxiliary field. 
In the following section we shall treat them as dynamical 
fields and generalize the ansatz of Eq.(9) to include the 
quanta of this field along with those of the fermions. We 
shall illustrate this for the symmetric case i.e. when there 
is no chemical potential mismatch for the two species of 
fermions. 



IV. DYNAMICAL CONDENSATE FIELDS AND 
GENERALIZATION OF THE BCS ANSATZ 



kinetic terms and the quartic interaction term, this La- 
grangian reduces to the one considered in section II. 

As before, we shall consider a state such that ($) = 
(j>o = (<&t) and investigate the fluctuations of the conden- 
sate field by defining the quantum fields $' = $ — 4>o, 
and $'t = $t _ O _ Then Cb reduces to 

- (V#' f )(V$') - A($'t$') 2 - 4A0 2 1 ($' t $') 

- Wo), (52) 

where, ~ means neglecting odd powers in <I>'s. Further, 
in the above, Vb(</>o) is the 'tree' level potential given as 



^o(0o) = K - tim + m 



(53) 



One can write down the corresponding Hamiltonian den- 
sities as 



with 



H = H f +H b + H hf (54) 



H t = J2 $ H a • V + /3m)tPi (55) 



As discussed earlier, we shall now treat the the scalar 
filed $ introduced in Eq.(3) as a dynamical field and 
rewrite the Lagrangian as 



where, 



C = C f +C b + C 



i>f 



0\„l,i 



Cf = -0 l (*7 M ^ - m + M7 jV' 



(48) 
(49) 



- n*, t *' t ) + *' t (-V 2 + mg)* / 

+ A($' t $') 2 +4A^($'t$') + T/ («/, ) 



(56) 



and TCbf = —£bf- Here, 11$/ (n$/t) is the conjugate 
momentum of the corresponding field $' ( ( f > '^). Similar 
to the fermion field operator expansion we can take the 
boson field operator expansions in terms of creation and 
annihilation operators. We have thus e.g. 



C b = (9o-^MB)$ t (9o+^MB)*-^$ t $ 

- (V$ t )(V$)-A($ t $) 2 (50) 



and, 



C bf = g\e^\ (VS 5 ^* + V>c7V$*) 



(51) 



We shall illustrate the effect of the dynamical bosonic 
field on the BCS-BEC crossover physics and hence will 
consider the case where there is no mismatch in the chem- 
ical potentials of the two condensing fcrmioionic species 
with a common chemical potential fi.. For the dynam- 
ical bosonic field now we have introduced the chemical 
potential [Ib which is equal to twice the fermionic chem- 
ical potential /i in equilibrium. We have also included 
a quartic term in the scalar field in Lb for the sake of 
completeness. In fact, this quartic term leads to non- 
perturbative corrections to the thermodynamic potential 
as we shall observe later. Clearly, in the absence of the 



$'(x,£ = 0) 
1 



dk- 



(2tt)3/2 J v /^(fc) 



= (a(k) + fot(_ k ))e* k - x (57) 



and 



n$/ (x, t = o) 



(2tt) 3 / 2 



dk 



w(jfc) 



(-fe(k)+a t (-k))e ik - x (58) 



Let us note that with the above expansion for 
the conjugate fields satisfying the quantum algebra 
[$'(x), 11$/ (y)] = id(x — y) leads to the usual commuta- 
tion relations for the creation and annihilation operators 
[a(k),at(k')] = S(k-k') = [&(k),&t(k')] for any arbitrary 
function w(k). 

With the operators for the scalar fields defined, we now 
gcncneralise the ansatz given in Eq.(9) and write down 
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the ansatz \Q(0,[i)b to include the effects of the boson 
field as, 



|Q(/3,M)>s=^C/ i W,M)>, 



(59) 



Here, similar to Eq.(9) and Eq.(8) the operator U B is 
given as, 



U B = cxp ( [ dkg(k)a?(k)tf(-k) - h. 



(60) 



and similar to Eq.(10) and Eq.(ll), the operator corre- 
sponding to thermal excitations of the bosonic fields IA B fx 
is given as 

Ufa = cxp(fdk(a'\k)a(-k)e a (k,fi) 

+ 6t(k)6(-k)0 6 (k,/*)) -h.c\ (61) 

This leads to e.g. 

(*'(x)*'(y)> 



l 



dk 



= ik-(x-y) 



(2tt) 3 J 2w(k) 
(cosh2g(k) + smh2g(k))(cosh 2 9 a + sinh 2 8 b ) 
= /(x-y,/3) 



(62) 



where, the expectation value is taken with respect to the 
state \Q,((3,h))b defined in Eq.(59). The extra functions 
g(k) as well as the thermal bosonic functions 9 a,b are to 
be determined as earlier by extremisation of the thermo- 
dynamic potential. The thermodynamic potential for the 
boson fermion system can be written as 



o, 



n + n B 



(63) 



Here, the fermionic contribution 51 has already been eval- 
uated in Eq.(23). The bosonic contribution £Ib is given 

by 



B = e B - [i-BpB - -5SB 



(64) 



The contribution of the first two terms is just the ex- 
pectation value of Eq.(56) with respect to the state in 
Eq.(60). The bosonic entropy density sb is given similar 
to their fermionic counterpart in Eq.(28), as [26], 



»B 



(27T) 



— rg y / dk( cosh 6 a In cosh 6 a 



sinh 9„ In sinh 6„ 



(65) 



Extremisation of the total thermodynamical potential 
with respect to the fermionic functions /(k), /i(k) and 
8* (k) leads to the same solutions for them as given in 



section II B. Extrcmising the bosonic function g(k) leads 
to the solution 



tanh2g(fc) = — 



M 2 



(66) 



with the quantity M 2 satisfying the temperature de- 
pendent 'mass gap' equation given as 



M z = m 
4A 



4A0 2 
dk 



(2tt) 3 J 2Vk 2 + M 2 
m 2 +4A(</> 2 + 1(0)). 



(cosh 2 8 a + sinh 2 #5 



(67) 



Here, 1(0) = 1(0,0), as given in Eq.(62). Similarly, min- 
imizing the thermodynamic potential with respect to the 
bosonic thermal functions yields 



sinh 9 a = 



n B (k) 



exp(E B - Ms) - 1 
for the boson particle distribution function and 



(68) 



sinh 9h 



1 



n B (k) (69) 



exp(E B +Ms) - 1 

for the boson anti-particle distribution function and 
Eb = Vk 2 + M 2 , with the temperature dependent 
mass M satisfying the self consistent mass gap equation 
Eq.(67). 

With all the functions in the ansatz state Eq. (59) now 
determined, the total thermodynamic potential can be 
written as, using Eq. (34) and Eq.(63), 



a 



flf + O s 



(70) 



where, the fermionic part of the thermodynamic potential 
is given by 



ty 



(2tt) 3 
2 

(2tt) 3 



(2e(k)-uj- -Q+)dk 



E 



dk 



In 1 



J-P^- } ) 



ln(l + e (-K s> ) 



(71) 



As compared to expression in Eq.(34), the above differs 
by the mass term m 2 ^ 2 , which is absorbed naturally in 
the bosonic part of the thermodynamic potential and the 
later is given as 



0, 



1 



(2tt) 3 
- 2\I 2 (0) 
+ V , 



dky/k 2 + M 2 



0(2n) 



2 

^s/^EMl-exp^)) 



(72) 
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where, the summation is over the bosons and antibosons 
with Ei = Eb — Hb, ^2 = Eb + Hb, and, Vq is the tree 
level potential 



V = (ml - nDcg + Ht 



(73) 



Finally, the extremisation of the total thermodynamic 
potential with respect to <j>o leads the superconducting 
gap equation 



M l 



Ms 



2A$ + V 



rfk 

(2^ 



cos 20_ 



uj- 



cos 26 , 
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— = -+4I 2 (A,Hs 
*R A 






X 



+ 4/ 1 (A, / u sc ) 



(75) 



(76) 



where, I\ and J2 are divergent integrals evaluated with 
a three momentum cutoff defined as 



- J_ /*' k|<A *_ 1- _^_ 
1 ~ (2^) 3 7 2|k| ~ A^o 8tt 2 



(77) 



and, 



|k|<A / , 

h = ^^r I dk 



m2c(2tt) 



2|k| 2V^L+1? 



! vr^i-i 



16tt 2 



/'■* 



(78) 



where, [i sc is the renormalization scale and A is the three 
momentum cutoff. Using the renormalization parameters 
the mass gap equation can be written in terms of finite 
quantities as 



M 2 =m 2 R + 4X R (cb 2 +If([3)) 



(79) 



with If(P) given as 



If (13) 



M 2 
16tt 2 
dk 1 



In 



M 2 



(2tt) 3 2E 



(n B (k) + n s (k)). (80) 



Similarly the bosonic part of the effective potential 
Eq.(72) in terms of the renormalized parameters, sub- 
tracting out the vacuum terms becomes finite and is given 
as 



(74) 

with, the mass M satisfying the mass gap equation 
Eq.(67) and cos 26^' 2 = (1-sin 2 6^(k)-sin 2 0^(k)) with with 
sin 0±(k) being the thermal distribution functions for 
the fermions defined in Eq.(31). This is the parallel of 
Eq.(35) where the condensate field was considered as an 
auxiliary field and there was no quartic coupling term for 
the scalar field. 

The expression for the bosonic part of the thermody- 
namic potential fl R in Eq.(72), is, however, affected by 
two types of divergences, one arising from the divergent 
integrals as vacuum terms (<po = at T = 0, /j, — 0) which 
are still to be subtracted, and, the other being the fact 
that the mass parameter M 2 is not well defined because 
of the infinities in the mass gap equation Eq.(67). This 
can be taken care of by defining the renormalized quar- 
tic coupling and the boson mass through the relations 
[21, 22] 



n B = v + v 1 + v 2 



(81) 



Vi = 



and, 



V = m R cb 2 + \ R <f>* - VbK + (A fl - A)$ ( 82 ) 



1 f 2 Af 4 ( M 2 



V 2 = -2X R I f {f3) 2 



(83) 

(84) 



The cut off dependence in the effective potential is still 
there in the last term of Eq. (82) which disappears in 
the limit of A — > 00 as the bare coupling A vanishes in 
that limit. In the present calculations, however, we keep 
the cutoff finite. We might note here that the difference 
between the bare and the renormalized quartic coupling 
can be written as 



4A K / 2 (A, Msc ) 

x - Xr -i-ax r i 2 (a^ sc ) Xr 



(85) 



where, I 2 is given in Eq.(78). Mostly in our numerical 
calculations we shall have the limit 4XrI 2 < < 1 in which 
case, the contribution of the last term in Eq.(82) is neg- 
ligible. 

Next, the fermionic part of the thermodynamic poten- 
tial Eq.(71) and the terms Vq, V 2 of the bosonic part can 
be combined using the gap equation Eq.(74) to yield a 
form similar to Eq. (41) as, 



n} = ii f + v + v 2 

2 



(2tt) 3 /3 



dk 



£--w- 



A 2 



£+ -UJ+ + 



2w_ t 2u> 

- 2 { In (l + e ("^- ; ) j + In (l + e (-H°)) } 

i 

- X R ^-(X R -X)4>t-2X R I f (p) 2 . (86) 

The number density equation, Eq.(37) now gets mod- 
ified due to the presence of dynamical condensate fields 
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having a chemical potential twice as that of the fermions 
and the number density equation becomes 



8/X0 2 , + 2 

1 



fsinh 2 6 a - sinh 2 6 b ) 



(2tt) 3 



dk 
— cos26»i 2 -— COS26* 1 ' 2 



dk (87) 



To discuss the crossover, we define the mass parameter 
mi as in Ref. [6] as 



2 2 a 2 

m i = m R - 4 3 



rfk 



(2tt) 3 Vk 2 + m 2 



and define the crossover parameter as 



A% 



(89) 



The gap equation Eq.(74) becomes 

m 2 -/4 = 2(\-\ R ) ( f> 2 -2\ R $~4\ R I f (f3) 



W 



dk 



cos 20 ' cos 26 , ' 



Vk 2 



(90) 



Let us note here that the superfluid gap equation has 
now contributions from the bosonic fluctuations through 
the contribution If(/3). For numerical calculations, we 
choose a given value of the bosonic mass parameter m R 
and solve the coupled number density equation and the 
superfluid gap equation Eq.(90) for the chemical poten- 
tial and the superfluid gap. At each stage of evaluation of 
the rhs of Eq.(87) and Eq.(90), the boson mass param- 
eter M 2 is calculated self-consistently solving Eq.(79). 
Throughout the numerical calculations we have chosen 
g = 2-\/2 and cutoff scale \j? sc = m R . This way we ob- 
tain the superfluid gap and the chemical potential for a 
given value of m R or equivalently, for a given value of x 
obtained through Eq.(88) and Eq.(89). These are then 
used to calculate any other thermodynamical quantities. 
In Fig. (10), we show the numerical solutions of such a cal- 
culation for the gap and the chemical potential for zero 
temperature as a function of the dimensionless order pa- 
rameter x/xq with xq = A 2 /(47r 2 ) for three values of the 
quartic coupling A_r=0, 0.5 and 2.0. A = 0, will corre- 
spond to the mean field results. As we might observe, 
the effect of the bosonic fluctuations are almost negligi- 
ble in the BCS regime (x/xq < 0) as well as near the 
unitary regime (x/xo ~ 0). The effects of the fluctuating 
field is seen at large values of the crossover parameter 
in the deep BEC regime, manifesting in a small reduc- 
tion of the superfluid gap. The magnitude of this reduc- 
tion increases with the quartic coupling. This is due to 




FIG. 10: Fermion chemical potential and gap in units offermi 
energy as a function of dimensionless crossover parameter 
x/xo . The solid (black), dotted (blue) and dashed (red) curves 
correspond to \r—0, 0.5 and 2 respectively. 



the following reason. It is clear from observing the gap 
equation Eq.(90) that, the effect of the positive quartic 
coupling leads to an increase in the bosonic chemical po- 
tential for the same value of the gap parameter. This 
leads to an increase of the corresponding value of the 
crossover parameter 'x' as may be clear from eq.(89) for 
the same superfluid gap parameter. This is also reflected 
in the number densities of the fermions and bosons which 
we show in Fig. (11). The change in chemical potential 
and the gap is appreciable from their corresponding mean 
field values only in the BEC regime. 

We then present the results for the critical temperature 
and the chemical potential as a function of the crossover 
parameter. This is done by setting A = = cf>o in the 
number density equation Eq.(87), in the gap equation 
Eq.(90) as well as in the mass gap equation Eq.(79) . 
The results are shown in Fig. 12. While the behaviour of 
the chemical potential is qualitatively similar to that at 
zero temperature, the critical temperature behaves sim- 
ilar to the gap at zero temperature with the crossover 
parameter. The correction to the critical temperature 
becomes significant in the BEC regime only and increases 
with the quartic coupling. In this case, again the reduc- 
tion of the critical temperature lies due to the increase 
in chemical potential due to thermal as well as vacuum 
fluctuations as may be seen in Eq.(90). With further in- 
crease of the quartic coupling, the critical temperature 
become less and less steep due to larger contributions of 
the thermal fluctuation of the the bosonic field. 

To see the effect of the fluctuations of the condensate 
field we also looked into the behavior of the gap as a func- 
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FIG. 11: Number densities of fermions and bosons in units 
of total number density. The solid (black), dotted (blue) and 
the dashed (red) curves correspond to Xr =0,0.5 and 2 respec- 
tively. 



tion of the temperature for different values of the quartic 
coupling parameter A_r. As the magnitude of the quartic 
coupling is increased, it was observed that the order pa- 
rameter A changes discontinuously at the critical temper- 
ature. A typical behavior for the solution of gap equation 
Eq.(90) is shown in Fig. 13 for Xr = 5. We note that near 
the critical temperature there are solutions to the gap 
equation but having larger values of the thermodynamic 
potential as compared to the normal matter. This is sug- 
gestive of a first order phase transition when the effect of 
fluctuations becomes large. One should however be care- 
ful in drawing conclusions from extrapolation to such a 
large value of the quartic coupling, because, although the 
result here is nonperturbative, it is limited by the ansatz 
for the ground state in Eq.(59). In this context we might 
remark here that gauge field fluctuations in color super- 
conductors changes the superconducting phase transition 
to a first order transition [23] . Similar observations were 
also made in Ref.[ll] in a fermion-boson model where the 
fluctuations were treated within a CJT formulation. 



SUMMARY 
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FIG. 12: Critical temperature T c and the chemical potential at 
T — T c in units of fermi energy as a function of dimension- 
less crossover parameter . The solid (black), dotted (blue) 
and dashed (red) curves correspond to Xr=0, 0.5 and 2 re- 
spectively. 
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We have considered here a variational approach to dis- 
cuss the ground state structure of system of two species 
of relativistic fermions with a mismatch in their Fermi 
momenta. An explicit construct for the ground state is 
considered to describe the two fermion condensates. The 
ansatz functions including the distribution functions are 
determined by an extremisation of the thermodynamic 



FIG. 13: Superfluid gap as a function of temperature for 
quartic coupling Xr=5. The dotted line corresponds unstable 
solutions which are to solutions of the gap equation but with 
a higher thermodynamic potential as compared to A = 0. 
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potential. 

The quadratically divergent gap equation is made log- 
arithmically divergent by subtracting out the vacuum 
contribution and the four Fermi coupling is related to 
the s-wave scattering length as in Ref.[7, 10]. Unlike the 
nonrelativistic case, the antiparticle degrees of freedom 
become important even for the case kf/m << 1, partic- 
ularly for large values of rj(= 1/fc/a). 

When the Fermi momenta of the two species are dif- 
ferent, we do not observe any gapless modes in the BCS 
regime. A nonzero chemical potential difference can sup- 
port a uniform BCS pairing with zero number density 
difference between the two species. Breached pairing so- 
lutions with two Fermi surfaces are also not observed. 
However, in the BEC region with p, < m, stable gapless 
modes are possible. The quasi particle can become gap- 
less for rj(= 1/kf-a) > 1.9 and even for larger value of 77, 
the quasi anti-fcrmions can also become gapless. Such 
gapless modes will be relevant for the transport coeffi- 
cients of the fermionic system. 

We have not calculated here the Meissner masses, or 
the number susceptibility to discuss the stability of differ- 
ent phases by ruling out regions in the parameter space of 
gap and chemical potential difference. Instead, we have 
solved the gap equation and the number density equation 
self consistently and have compared the thermodynamic 
potentials. In certain regions of the chemical potential 
difference and the coupling, we have multiple solutions 
for the gap equation. In such cases, we have taken the 
solution which has the lowest thermodynamic potential, 
ensuring that it corresponds to a minimum. In the deep 
BEC region, the phase transition from BCS to gapless 
phase is a second order phase transition with the order 
parameter decreasing continuously, while the transition 
from the gapless phase to the normal matter phase is a 
first order transition as the difference in the densities of 
the two condensing species is increased. 

The results obtained here are of course limited by the 
choice of the ansatz. We have not considered here other 
nonuniform ansatz leading to the LOFF phase or the 
crystalline phase [14, 15, 32, 33]. The results obtained 
here might nevertheless be regarded as a reference so- 
lution with which other numerical or analytical results 
obtained from more involved ansatz for the ground state 
may be compared. Though it is not apriori obvious, the 
results, obtained through the simple variational ansatz 
for the ground state, regarding the phase structure for 
this purely fermionic theory turn out to be similar to 



those of a Bose fermi model of Rcf. [6] treated within a 
mean field approximation. 

We have also considered the effects of the fluctuations 
by treating the condensate field as a dynamical bosonic 
field in a model with quartic self interactions of the boson 
field. The BCS ansatz was modified to include the quanta 
of the fluctuating field along with the usual fermion pairs. 
Such an ansatz gave rise to a supcrfluid gap equation that 
includes the effect of condensate fluctuations. In the eval- 
uation of the superconducting gap the scalar field mass 
gap was also calculated self consistently. This leads to a 
decrease of the critical temperature in the BEC regime. 
We also observed that the superfluid transition could be 
first order for larger quartic coupling with the effect of 
the condensate fluctuations becoming larger. The present 
ansatz for the ground state leads to the result arising from 
a summation of an infinite series of bubble diagrams for 
the scalar field. However, this does not include the effect 
of 'sunset' type diagram. Inclusion of such diagrams has 
been successfully done recently within a C JT formulation 
in Rcf. [11]. 

We can generalize this toy model to a more realistic 
model like Nambu Jona-Lasinio model, to describe the 
possibility of relativistic BCS-BEC crossover in quark 
matter phase diagram. Further, the effect of charge neu- 
trality conditions, as appropriate for matter in the core 
of neutron stars, on relativistic BCS-BEC crossover will 
be important. Some of these calculations are in progress 
and will be reported elsewhere. 
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